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Abstract 



A new approach is presented for the problem of optimal impulsive rendezvous of a 
spacecraft in an inertial frame near a circular orbit in a Newtonian gravitational field. 
The total characteristic velocity to be minimized is replaced by a related characteristic- 
value function and this related optimization problem can be solved in closed form. The 
solution of this problem is shown to approach the solution of the original problem in 
the limit as the boundary conditions approach those of a circular orbit. Using a form 
of primer-vector theory the problem is formulated in a way that leads to relatively 
easy calculation of the optimal velocity increments. A certain vector that can easily be 
calculated from the boundary conditions determines the number of impulses required 
for solution of the optimization problem and also is useful in the computation of these 
velocity increments. 

Necessary and sufficient conditions for boundary conditions to require exactly three 
nonsingular non-degenerate impulses for solution of the related optimal rendezvous 
problem, and a means of calculating these velocity increments are presented. If neces- 
sary these velocity increments could be calculated from a hand calculator containing 
trigonometric functions. A simple example of a three-impulse rendezvous problem is 
solved and the resulting trajectory is depicted. 

Optimal non-degenerate nonsingular two-impulse rendezvous for the related prob- 
lem is found to consist of four categories of solutions depending on the four ways the 
primer vector locus intersects the unit circle. Necessary and sufficient conditions for 
each category of solutions are presented. The region of the boundary values that ad- 
mit each category of solutions of the related problem are found, and in each case a 
closed- form solution of the optimal velocity increments is presented. Some examples 
are simulated. Similar results are presented for the simpler optimal rendezvous that re- 
quire only one-impulse. For brevity degenerate and singular solutions are not discussed 
in detail, but should be presented in a following study 

Although this approach is thought to provide simpler computations than existing 
methods, its main contribution may be in establishing a new approach to the more 
general problem. 



1 Introduction 



It has been a convenient practice to model spacecraft orbital and trajectory problems having 
relatively high thrusts over short intervals of time by discontinuous jumps in velocity, but 
retaining continuity of the position vector at the time at which the discontinuity in the 
velocity appears. Problems that are modeled in this way are called impulsive orbital or 
impulsive trajectory problems. The impulsive problems considered here are based on 
the restricted two-body problem, (i.e. a particle of mass in a Newtonian gravitational field 
emanating from a point source). These minimization problems consist of the determination 
of a finite set of velocity increments Av 1? ... Av^ and the values of the true anomaly 0\ ... Ok 
at which they are applied, to minimize the total characteristic velocity E|Avj| subject to the 
two-body equations of motion, the initial position t(9q) — r o, initial velocity v(#o) = Vo, the 
terminal position r(0f) = ry and the terminal velocity v(#j) = vy where Oq an d Of denote 
respectively the initial and terminal values of the true anomaly. If r , v , ry, vj are specified 
points in M 3 the minimization problem will be called an optimal rendezvous problem. 
If r , v or Tf, Vf or both are arbitrary on a nontrivial arc of a fixed two-body orbit, the 
minimization problem will be called an optimal transfer problem. Clearly, any solution 
of an optimal transfer problem found at the points r(0 ), v(0 o ), r (0f), v(0/) also defines a 
solution of an optimal rendezvous problem having those end conditions. 

Apparently the first significant publication of an impulsive orbital problem was the 
Hohmann Transfer [1], an elliptical arc, tangent to two circular orbits, that appeared in 
1925. This was followed in 1929 by the Oberth Transfer [2], a two-impulse transfer connect- 
ing a circular orbit to a hyperbolic orbit through an ellipse with apsides touching the circle 
and the center of attraction. The second impulse at (or near) the center of attraction sends 
the craft into hyperbolic speed. 

Some of the earliest studies of orbital maneuvers were done by Contensou [3, 4] and 
Lawden [5]. Surveys of much early work were done by Edelbaum [6], Bell [7], Robinson [8], 
and Gobetz and Doll [9]. An excellent introduction to the subject is found in the book by 
Marec [10]. 

In the fifties it was discovered that the total characteristic velocity can be reduced over 
that of the Hohmann Transfer by the addition of a third impulse if one of the circular 
orbits is much larger than the other [11, 12, 13]. This might lead one to expect that some 
minimum total characteristic velocity problems for a single Newtonian gravitational source 
would require at least three impulses. In 1965 a paper appeared by Marchal [14] showing 
that a bounded version of the bi-elliptic transfer was indeed a three-impulse optimal transfer. 

Recently this problem has been revisited by Pontani [15] who through calculus and a 
simple graphical technique separated the optimal two-impulse solutions (Hohmann) from 
the optimal three-impulse solutions (bi-elliptic). 

Although three-impulse solutions to the optimal rendezvous problems and the optimal 
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transfer problems exist, they are somewhat rare in the literature. Ting [16] in 1960 and 
Marchal [14] in 1965 discovered that at most three impulses are sufficient to solve the optimal 
rendezvous or transfer problem. There are solutions having more than three impulses, but 
these solutions are either degenerate or singular solutions, and also have a three-impulse 
solution containing the same total characteristic velocity as those having more than three- 
impulses. 

Apparently the required number of impulses for optimality is dependent on the way 
the equations of motion are modeled. In relative-motion studies in which the equations of 
motion are linearized about a point in Keplerian orbit as many as four impulses can be 
required for optimality of the planar problem. It has been shown by Neustadt [17], and 
others [18-20] that the maximum number of impulses needed for linear equations is the 
same as the number of state variables, that is, four for planar problems. In 1969 Prussing 
[21] displayed optimal four-impulse rendezvous maneuvers using equations linearized about 
a point in circular orbit. This problem was revisited much later by Carter and Alvarez 
[22]. Calculations of several optimal four-impulse rendezvous maneuvers and a degenerate 
five-impulse rendezvous maneuver were presented by Carter and Brient [20] in 1995 using 
equations linearized about elliptical orbits. 

Optimal velocity increments in impulsive problems are generally calculated in one of 
three ways: The solution of a Lambert's problem, primer- vector analysis, or parameter 
optimization such as through a nonlinear programming algorithm. It is common practice to 
combine these methods in the solution of specific problems. In the present paper we employ 
a form of primer-vector analysis. 

The primer vector was introduced by Lawden [5] as the part of the adjoint vector cor- 
responding to the velocity that satisfied certain necessary conditions for optimality of a 
trajectory. Lion and Handelsman [23] summarized these necessary conditions and showed 
how to construct improvement in the velocity increments based on primer vector analysis. 
Prussing [19] showed that the necessary conditions of Lawden and Lion and Handelsman are 
also sufficient if the equations of motion are linear. 

The present paper uses equations of motion that describe the restricted two-body prob- 
lem, but the analysis is confined to near-circular orbits. This limitation is caused by an 
approximation which is highly accurate when the radial speed is small compared with the 
tangential speed. This simplification has been used successfully by the authors in studies 
that involve atmospheric drag [24, 25]. We replace the total characteristic velocity which is 
the original cost function by a related cost function which is useful near a nominal circular 
orbit, and present closed- form solutions to this related optimization problem. These results 
should be more accurate than those based on the approximate Clohessy- Wiltshire equations 
[28] because the linearization used herein is not approximate. We remark that the trajec- 
tories described between impulses are not approximations, and the solutions we present are 
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optimum in terms of the model that we use. We may therefore use the terminology "optimal 
solutions". These velocity increments, however, are not accurate unless the trajectory is 
near a circular orbit. 

This paper uses a form of primer- vector analysis in a novel way. It has been well known for 
ages that in the restricted two-body problem, the equation of motion involving the magnitude 
r of the position vector r can be transformed to an equation linear in its reciprocal. This 
forms the basis of a set of linear equations that describe the motion of a particle of mass 
in the transformed variables. The impulsive minimization problem is then formulated in 
terms of the resulting linear equations, and is amenable to a specific theory of impulsive 
linear rendezvous developed by Carter and Brient [26, 20, 27]. The definition of the primer 
vector for linear equations is in a more general form than that of Lawden[5] or Lion and 
Handelsman[23], but the more general form is needed for use with the necessary and sufficient 
conditions for solution of the impulsive minimization problem studied in this paper. 

In this new formulation the state variables r, the radial distance from the center of 
attraction, f its time rate of change, and , the true anomaly have been replaced by variables 
Vii ViiVz to be defined in the next section, called transformed variables. The transformed 
state vector y = (i/i, 2/2, 2/3) T satisfies a linear differential equation. It is straightforward to 
define a transformed state transition matrix. The primer vectors are shown to be ellipses for 
this problem. It will be shown that for three-impulse solutions and two-impulse solutions the 
primer vectors are completely determined by the initial true anomaly and the terminal 
true anomaly Of. 

Given 0q and Of the initial transformed state vector yo and the terminal transformed 
state vector y/ are used to define a generalized boundary point z/(yo,yf). The geometric 
structure of the set of generalized boundary points associated with an optimal primer vector 
is a convex cone [20]. Given, 9 and Of these convex cones partition the set of generalized 
boundary points into a simplex of convex conical sets. The three-dimensional cones contain 
the generalized boundary points that admit three-impulse solutions. Points on the boundary 
of such a cone admit degenerate three-impulse solutions, that is, one or more of the three 
velocity increments is zero. Although three-impulse solutions are somewhat rare in the 
literature for the restricted two-body problem, this analysis shows that there are plenty of 
them near circular orbits. It follows also that the two-dimensional cones define areas that 
admit two-impulse solutions, and the one-dimensional cones, that is, the straight lines admit 
the one-impulse solutions. Of course the vertex of any cone, that is, the origin, admits only 
the zero-impulse solution, a coasting trajectory. 

For the original problem having boundary conditions in the vicinity of a nominal circular 
orbit, near optimal velocity increments can be obtained in closed form. For more general 
boundary conditions the velocity increments calculated from the related problem could be 
useful as good initial guesses in a numerical optimization program, and should obviate the 



5 



need to solve Lambert's problem. 

The paper begins with the general formulation of the problem and the related problem 
for boundary conditions near circular orbits, then presents some analysis and simulations 
of non-degenerate three-impulse solutions and various kinds of non-degenerate two-impuse 
solutions of the related problem and precisely determines the sets of boundary conditions 
corresponding to each type,followed by some results on one-impulse solutions. For brevity 
detailed discussion of degenerate and singular solutions is postponed to a later work. 



2 Building the Model 

The equation of motion of a particle in a Newtonian gravitational field about a homogeneous 
spherical planet is 

r = -4r (2.1) 

where r is the position vector with respect to the center of the planet, \i is the product 
of the universal gravitational constant and the mass of the spherical planet, the upper dot 
indicates differentiation with respect to time t and r = |r| = (r ■ r) 1 / 2 where the inside dot 
denotes the usual inner product, and |r| represents the Euclidean norm of the vector r. 



2.1 The Optimal Impulsive Rendezvous Problem 

In polar coordinates (12. li becomes 

rr + 2f6 = (2.2) 

r- r 6 2 = -^ (2.3) 

where 9 represents the true anomaly. The radial and transverse velocities are denoted re- 
spectively by v r = r and vg = r9. Throughout we shall use subscript notation to refer to 
any variable x at time as = x(ti). 

We shall consider the addition of velocity increments Avj = (Av ri , AygJ 7 at time t{ for 
i — 1, k such that 

lim r(t) = lim r(t) + Av 4 . (2.4) 

t— >u+ t—>ti— 

In this formulation r and 9 are continuous everywhere and the differential equations (12.21) 
and (12.31) are satisfied everywhere on an interval to < t < tf except at the instants where 
the velocity satisfies the jump discontinuities (12.41) . At the ends of the time interval, (12.41) 
becomes 

lim r(t) = r(t ) + Avi, t x = t (2.5) 

r(t/)= lim r(t) + Av fc , = (2.6) 

t-*tf- 
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The rendezvous problem requires that the initial conditions 



r(t ) = r , 8(t ) = 6 t 



o 



(2.7) 



r(t ) = Vro , v e (io) = v flo 

and the terminal conditions, 



(2.8) 



r(t f ) = v f , 6(tf) = e f 
*(*/) = v r/, v*(t/) = v« 



(2.9) 
(2.10) 



be satisfied. The total characteristics velocity is defined by 



fc 



c = 



2>v, 



(2.11) 



The optimal impulsive rendezvous problem is stated as follows: 

Find a positive integer k, a finite set {ii, ...t&} on the interval {to < t < and a 
set of velocity increments {Avi, ... Av^} such that the differential equations ( 12. 2ft and f ]2. 3f) 
are satisfied except on the finite set {t\, ...tfc} where ( 12.41) - ( 12761) are satisfied, the boundary 
conditions (l2.7p -( l2"TTUj) are satisfied, and the total characteristic velocity ( 12. lip is minimized. 

Similarly one can state an optimal impulsive transfer problem by replacing the require- 
ment that the initial conditions, as described by ( 12. 7p . (12. 8p and the terminal conditions as 
described by (12.91) . ( 12.101) be fixed with the requirement that either or both be contained on a 
Keplerian orbit. Mathematically, the optimal impulsive transfer problem is not significantly 
different from the optimal impulsive rendezvous problem. In fact, an optimal solution sat- 
isfying the fixed end conditions ( 12 . 7|) - fl 2 . 1 [) also defines an optimal solution of the transfer 
problem from the Keplerian orbit containing the fixed point described by (12.71) . (I2.8P to the 
Keplerian orbit that contains the fixed point described by (12.91) . (12.101) if the time interval 
is sufficiently large. 



We sketch the derivation of the orbit equation which has been well known for many years. 
Multiplying (12. 2p by r, we have the derivative of 



The constant h is the angular momentum. We may use (12.121) and the chain rule to change 
the independent variable in (12.31) to 9. We use primes to indicate differentiation with respect 
to 9. By this change of variable and noting from (12.121) that 9 is monotone in t, the fourth 



2.2 The Orbit Equation 



r 2 9 = h. 



(2.12) 
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order system (I2.2p . (12.31) on the interval to < t < tf is replaced by the following third order 
system on the interval 9q < 9 < 9 f\ 

r(9)r"(9) - 2r'(9) 2 = r(9) 2 - fir(9) 3 /h 2 (2.13) 

h'{9) = 0. (2.14) 

2.3 Transformation to Linear Equations 

Employing the well-known change-of- variable y — 1/r in the orbit equation ( 12 .131) and ma- 
nipulating some symbols, we obtain the linear differential equation: 

y" + y = ^/h 2 . (2.15) 

We now assign the state variables: 

2/1 = 2/ (2.16) 

V2 = y' = -r/h (2.17) 

y 3 = fi/h 2 . (2.18) 

In terms of these state variables, the equations (12.131) . (12.141) are transformed into the fol- 
lowing set of linear equations: 

y[ = y 2 (2.19) 

2/2 = -2/i + 2/3 (2.20) 

V'z = 0. (2.21) 

The initial conditions are 

Vi(0 o )=y lo =- (2.22) 

2/ 2 (0o) = 2/20 = ~m)/ho = -— (2.23) 

r v eo 

M) = «- = & = (2 ' 24) 

and the terminal conditions are 

!&(*/)= 01/ = -• (2-25) 

= 2/2/ = -r(*/)/*/ = -^f- (2.26) 
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where the subscript refers to the initial value of a variable and the subscript / refers to 
the terminal value. Formulas for recovery of the original variables are 



l/2/i (2.28) 



e = yl^fm (2.29) 

r = -yi^lVz (2-30) 
vo = ViVfi/ya- (2-31) 

2.4 Restatement of the Optimal Impulsive Rendezvous Problem 

At 9i, i — 1, k we increment the velocity r(9i) by Avi = (Av n , AvgJ 7 so that (I2.4p -f l2l6|) 
are satisfied. This causes an instantaneous jump in the state variables y 2 and y 3 but y 1 is 
continuous. The increments in y 2 and y 3 caused by the increments in the velocity at 9i are 
determined respectively from ( 12.17ft and ( 12.181) . 

Ay 2i = —r(Av ri - ^Av 6i ), i = l,...,k (2.32) 

hi v Bi 

Ay 3l = -^Av 9i , i = l,...,k. (2.33) 

We now consider the imposition of a velocity increment Avj = (Av ri , Av0 t ) T at 9i for 
i = I, ...,k on the interval 9q < 9 < 9f such that 

lim y 2 (9) = lim y 2 (9) + Ay 2u lim y 3 (9) = lim y 3 (0) + A^. (2.34) 

If 9i is an end point (I2.34p becomes 

lim y 2 {9) = y 2 {9 ) + Ay 2i , lim y 3 (9) = y 3 {9 ) + Ay 3l (2.35) 

0->-e + 9-+8 + 

y 2 (9 f ) = lim y 2 (9) + Ay 2l , y 3 (9 f ) = lim y 3 (9) = Ay 3i . (2.36) 

u — >u j — u — >6j- — 

The optimal impulsive rendezvous problem can now be restated as that of finding a 
positive integer k, a finite set K = {9i,...,9k} on the interval 9$ < 9 < 9f, and a set 
of velocity increments {Av^ G M?,i = l,..,k} to minimize the total characteristic velocity 
( 12.1 ip subject to the linear differential equations (I2.19p - (l2.2ip that are valid except on the set 
K where f !2.32[) - (I2.36j) are satisfied, and that satisfy the boundary conditions ( I2.22p -( |2~2~T|) . 
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2.5 Transformed Velocity Increments and the Related Problem 

We transform the velocity increments Av Ti and Av$ i , i = 1, k as follows: 

AV U = 1 ( Av n - ^AvgA (2.37) 



AV 2l = ^Av 6i . (2.38) 

With this transformation, we define the increments Ayj G M 3 , i = 1, k from (I2.32p and 
(I2.33p . caused by the velocity increments as 

A yi = 5AVj (2.39) 

where AV { = (AVu, AV 2i ) T , i = 1, k and 

( °\ 

B = -1 . (2.40) 
V -2 J 

For i = 1, k the expressions ( 12.34j) -( 0.36[) . in vector notation become 

lim y(6) = lim y(0) + Ay,, 9 < 9 { < 9 f , (2.41) 

— >6 t -\- 6 — >6-i — 

lim y(e)=y(e ) + A yi , 9, = 9 , (2.42) 
y(9 f )= lim y(0) + A yi , 0< = 0,. (2.43) 
In vector form the boundary conditions ( 12.22l) -( |2.27l) are written 

y(0 o ) = y (2.44) 
y(0/) = y/ (2-45) 

where y = (yio, 2/20, Z/3o) T and y/ = (y lf , y 2f , y 3 f) T are defined from the respective right- 
hand sides of (EQ2D - ([01 and (1235) - (157711 . For £ AT the linear equations (l2~T3D - ([2T2Tj) 

are 

y = ^4y (2.46) 

where 

1 

A= | -1 1 | . (2.47) 
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We now replace the cost function (12. lip by the related cost function 



k 

C = ^2\AVi\. (2.48) 

i=l 

The related optimal impulsive rendezvous problem is that of finding a set K = {9i, ...,6k} 
on the interval 9 < 9 < 9f and velocity increments AVi, AV& G M. 2 that minimize 
the related cost function C subject to the linear differential equation (I2.46P that is satisfied 
everywhere on the interval except the set K where the velocity increments AVi, AV& are 
applied subject to ( 12.39[) and (12. 4 II) - ( 12^31) . and the initial condition (12.441) and the terminal 
condition (I2.45j) are satisfied. 

It has been shown [20] that a solution of this problem exists under appropriate control- 
lability conditions. 

Since the differential equation (12.461) is linear, it is known [17-20] that it is sufficient to 
set k = 3. For this reason we set k = 3 for the remainder of this paper. This agrees with 
the results of Ting [16] and Marchal [14] although the linearized model used by Prussing [21] 
allows k — 4. 



2.6 Necessary and Sufficient Conditions 

A fundamental matrix solution $ associated with a matrix A satisfies the matricial differen- 
tial equation $' = A$. It is not difficult to see that 

cos 9 sin 9 1 

<l>(0) = j -sinfl cos# J (2.49) 
1 

is a fundamental matrix solution associated with (I2.47p . The inverse of $ is 

(cos 9 — sin 9 — cos 9 \ 
sin# cos# -sin9 . (2.50) 
1 J 

Lawden's definition [5] of the primer vector is not adequate for the work herein. We use 
the definition [20,26]: 

p(0) = R(9) T X (2.51) 

where A 6 R 3 and 

R{9) = <&-\e)B. (2.52) 
We briefly review some previous results [20]. 
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If AV, is a velocity impulse at 0< it follows from f l2~39]) . (I2T46]) . and f lCT|) -f l2T43|) that 

y (0) = S^)^)- 1 ^ + 5AV,) (2.53) 

on the interval 8i < 6 < 9 i+ i if there is a velocity impulse at 0j+i> otherwise 9^ < 9 < Of. If 
there are a total of k impulses it follows from (12T39D . f l2T4TD -f l2^4l ) and (12331) that 

y(0) = $(0)$(0 o )- 1 yo + $(^) ^ $(^) _1 5AV,, (2.54) 

i=i 

where 9^<9<9j. Applying the terminal condition ( I2.45P we obtain 

k 

y f = $(^)$(0o)" 1 yo + H0 f ) my'BAV,. (2.55) 

i=l 

Introducing the vector z/, we write this expression as 

k 

z f = J2®(8i)~ lBA V i (2.56) 

i=i 

where all of the information about the boundary conditions is contained in the definition: 

z f = ^(9 f )- 1 y f -^9 )- 1 y . (2.57) 

Setting aii = |AVi|, i = 1 ... k we have all of the terminology needed to state the following 
result from previous work. [20] 

Theorem: For a minimizing fc-impulse solution of the related optimal impulsive ren- 
dezvous problem, it is necessary and sufficient that 



AVi = 0, 


or AVj = -p(9i)ai, i = l,...,k 


(2.58) 




oti > 0, i = 1, . . . , k 


(2.59) 




k 

y *r t R(6 i )p(6 i )a i = -z f 
i=i 


(2.60) 


AV,= 


= 0, or \p(6i)\ = 1, i = l,...,k 


(2.61) 


= 9 , or 


|p(0i)r = O, or 9i = 9 f , i = l,...,k 


(2.62) 




\p(0)\<h <9<9 f . 


(2.63) 



Proof: See Sec 3.2 and 3.3 of Reference 20. 

In this theorem k is the specified number of impulses and can be any non-negative integer 
sufficiently large that fl2.58p ~ fj2.63p are satisfied. It has been shown [17-20] to be unnecessary 
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for k to be larger than the dimension of the state space (i.e. dim zj which is 3 for the present 
problem), although there can be degenerate solutions where k > dimzf. A degenerate 
solution is a minimizing solution for which there is an equivalent minimizing solution where 
some on = in (12. 59ft . If a solution is degenerate an equivalent minimizing solution exists 
having fewer than k impulses so that, effectively, the number k could be reduced. 

For many boundary conditions there can be minimizing solutions for k < dimzf. For 
this reason we shall set k equal to the number of elements in the set 

# = {0|| P (0)| = o, e <e<e f } 

if this set is finite. Since the function f{9) = p T {9)p(9) — 1 is analytic either K is finite (i.e. 
K = . . . , 8k}) or else K is the entire interval 9q < 9 < Of. In the latter case a minimizing 
solution is called singular and k and 0\, ... ,9k are arbitrary as long as k is large enough 
that (|2.58p - (l2.63p are satisfied, as they will be if k = dimzf. 

In the present paper we shall isolate those boundary conditions where k = 3 and inves- 
tigate minimizing non-degenerate three-impulse solutions. In the following paper we shall 
investigate the non-degenerate two-impulse and one-impulse solutions. In the final paper we 
shall present some results on singular and degenerate solutions. 

We shall introduce some terminology based on (12 .62 p . We refer to 9 as an initial value 
of 9 and 9f as a terminal value of 9. If an optimal impulsive rendezvous has a velocity 
increment assigned at 9q (i-e. AVj = AVo) it is called an initial velocity impulse; if it has a 
velocity increment assigned at 9f (i.e. AVj = AVj) it is called a terminal velocity impulse. 
A value 9\ is called a stationary value if = and p(9i) = 1. A velocity increment 

AVj assigned at a stationary value 9i is called a stationary velocity impulse. If the one-sided 
derivative of \p{9)\ is zero at 9 , then 9 is both an initial value and a stationary value, and 
AV is both an initial and a stationary velocity impulse. Similarly, if the one-sided derivative 
of \p{9)\ is zero at 9f, then 9f is both a terminal value and a stationary value and AVy is 
both a terminal and a stationary velocity increment. 

2.7 Optimal Impulsive Rendezvous near a Nominal Orbit of Low 
Eccentricity 

Given 9q, yo, 9f, y/ an analytical solution will be presented for the related optimal impulsive 
rendezvous problem based on the necessary and sufficient conditions recently stated. It is 
possible to display the velocity increments, their points of application, and the actual arcs of 
Keplerian orbits that minimize the related cost function subject to (I2.39p -( |2.47p . A practical 
limitation of this result is that it is the total characteristic velocity (12. lip that we seek to 
minimize, not the related cost function (12.481) . 
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We shall show that the total characteristic velocity is approximately proportional to 
the related cost function and approaches exactness as y and y/ approach end conditions 
associated with a circular orbit. 

Consider a state vector y(9) associated with a nominal Keplerian orbit where y(#o) = yo 
and y(0f) = yj. Since there are no impulses in the nominal orbit, (I2.55f) shows that 

y/ = $ (fl/Wo) -1 ^- (2.64) 

Comparing an impulsive orbit with a nominal orbit, 

k 

y(0) - y(0) = W0o) _1 (yo - yo) + mr'BAV, (2.65) 

i=l 

where k = 1,2, or 3. From ( I2.56P and ( I2.57P we observe that 

z f -z f = $(^ / )- 1 (y / - y f ) - H9 )-\yo - yo) (2.66) 

where z/ = from (I2.56p . This shows that zj can be made arbitrarily small by making yf 
sufficiently near ff and y sufficiently near y . Upon multiplying both sides of f !2.56|) by a 
very small positive number, we see that we can select e > and |AVj| < e/3, i = 1,2,3 
where e is made arbitrarily small by selecting Zf sufficiently small. In this manner we can 
make the cost C < e and the difference (I2.65P in the trajectories arbitrary small. 

We shall choose the nominal trajectory y(6) to be a circular orbit or a low-eccentricity 
elliptical orbit. Picking end conditions sufficiently near the end conditions of the nominal 
orbit we establish a sufficiently small bound on AVj, % = 1,2,3 (i.e. there is a sufficiently 
small positive number M such that |AVj| < M for i = 1, 2, 3). This establishes a bound N 
on |At>#. |, for % = 1, 2, 3 (i.e. there is a number N > such that lAt^J < iV for i = 1,2, 3). 
We have shown that the difference (12.651) from the nominal circular or low-eccentricity orbit 
can be made as small as possible. This and (12.171) show that \v Ti \ is much less than [vg^ near 
a low-eccentricity nominal orbit. The conclusion of these arguments is that the second term 
in the parenthesis on the right-hand side of (I2.37P can be made arbitrarily small by picking 
the end conditions sufficiently near those of the nominal circular orbit. 

One way to select a radius defining a nominal circular orbit is by setting f = (ro +r/)/2. 
It is not difficult to show that the angular momentum h of a circular orbit of radius f is 

h = (fif) 1/2 . (2.67) 

Using k = 3 we now show that the total characteristic velocity (12.111) can be approximated 
to any desired accuracy using the related cost (I2.48P by picking the end conditions sufficiently 
near these of a nominal circular orbit. 
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It follows from f[2T3Tj) that for i = 1, 2, 3 
lAv; - hiAVA = 



v 



' - 2 I f" ; r-' 

+ 1 /L? J 



lV2 



(2.68) 



Given any e > 0, , there are end points sufficiently near end points of a nominal circular 

(2.69) 



orbit such that for i = 1,2,3 



< 



6N' 



fir jATj 



< 



6N 



It follows from ( 12 . 6 7|) and the fact that |Au^| < N for i = 1, 2, 3 that 

|Av< - hi AV,| < e -. 

6 



(2.70) 



To complete the arguments, we note that 

3 3 3 3 3 

|^|Avi|-^^|AVi|| < J2 lAvi-MVil < ^|A Vi -/iiAVi|+^|/ii-/i||AVi|. (2.71) 



8=1 



i=l 



i=l 



8=1 



8=1 



Picking initial and final conditions sufficiently near the nominal circular orbit, for each 
i = 1,2,3 we have 



\hi — h\< 



6M 



This establishes the final result: 



^lAvil-Ji^lAVi 



< e 



(2.72) 



(2.73) 



8=1 



8=1 



For trajectories near a nominal orbit of low eccentricity we have the approximation 



^|Av,| = /i]T|AV,,|. 



(2.74) 



8=1 



8=1 



Since h is a constant we may use the cost 



(2.75) 



8=1 



in these situations. After solution of this related problem one can calculate the actual 
impulses Avj in terms of AVj through (I2.37P and (I2.38P where v n and v$ i are recovered 
from (I2.28P - (I2.3ip . If the end conditions are very close to a nominal circular orbit one can 
simplify and approximate (I2.37P and f |2.38|) by 



Av, = h&Vi, i = 1, 2, 3. 



(2.76) 
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3 Primer Vector Analysis 



The primer vector can be a very useful tool [5, 23, 19, 20] in the analysis of optimal impulsive 
problems. In this section we survey the various geometric arrangements of primer vector loci 
for the problem addressed, its degeneration in the case of singular solutions, and its use in 
the calculation of three-impulse trajectories. 



3.1 Geometry of Primer Vector Loci 

From (I2.40p . ( 12.50jl . (12. 52ft and (12.5ip the primer vector is determined in terms of the vector 
A. Writing A T = (Ai,A 2 ,A 3 ) and p(9) T = \p 1 {9),p 2 {9)) it follows that 

Pi(0) = Ai sin9- A 2 cos 9 (3.1) 

p 2 (0) = 2Ai cos + 2A 2 sin^- 2A 3 (3.2) 

We set Ai = Acos0 and A 2 = Asin</> where A = (\\ + X 2 Y^ 2 - The primer vector is therefore 
described by 

Pi(0) = Asin(0-0) (3.3) 

p 2 (6) = 2A cos(0 - 0) - 2A 3 (3.4) 

If A 7^ then the primer vector loci are ellipses. The equations (13. 3p and (13.41) can be 
combined and put in the standard form: 



pj (pa + 2A 3 ) 2 _ 

A2 + (2A) 2 " ^ 

The major axis is 2A, the minor axis is A, and the center is at (— 2A 3 ,0) where the p 2 - 
axis is the ordinate and the pi-axis is the abscissa. According to (I2.6ip it is necessary that 
optimal velocity impulses are applied at the values of 0, described by the intersections of 
the ellipse (13. 5 p and the unit circle. Also (I2.63P shows that the primer vector arc must be 
contained inside the unit circle. This is depicted graphically in Figures 1-8 for various types 
of three- impulse, two-impulse, and one-impulse solutions. The arrows indicate the points of 
application of the velocity increments. 

If 9 is any value of the true anomaly we let 9 = 9 — <fi. Similarly for a subscript i we let 
9i = 9i- (p. We therefore write (Q and (E3D as 

Pi(0) = Asin0 (3.6) 

p 2 (Q) = 2A(cos0- k) (3.7) 

where 

k=^. (3.8) 
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If ( 13. 6 p and ( 13. 7p are associated with an optimal solution for a generalized boundary point 
Zf and optimal velocity increments are AVi, % = 1, k, we see that 

Pl (e) = -Xsm6 (3.9) 

p 2 (0) = -2A(cos#- k) (3.10) 

are associated with the generalized boundary point — Zf and optimal velocity increments are 
—AVi, i — 1) •••) fc. This latter solution for the generalized boundary point — Zf will be called 
the antipodal solution relative to the original solution associated with Zf. Obviously if the 
original solution has been determined, the antipodal solution is available without calculation. 
Figure 2 depicts the locus of a primer vector that is antipodal to that presented in Figure 1. 
This concept is not found much in the literature although Breakwell does refer to "mirror 
image" trajectories [29]. 



3.2 Singular Solutions 

An arc of a trajectory is called singular if \p(9) \ = 1 identically on an interval. If A = 0, then 
( 13~3|) . ((331), and (l23T|) show that A 3 = ±1/2 and 

Pi(0) = O (3.11) 

p 2 {6) = ±1 (3.12) 

everywhere on the interval 6$ < 6 < Of. There are no other values of Ai, A2, or A3 for which 
singular solutions occur. The two constant primer vectors described by (13. lip and (13 . 121) 
will define two optimal singular solutions on Oq < < Of. From ( I2.58p - (l2.60p we find that 

AV, = ±( ° J , i = l,...,k (3.13) 

k 

^oacosBi = ±z h /2 (3.14) 

i=l 
k 

Qj sin Oi = =FZ/ 2 /2 (3.15) 



i=l 

k 



J2<Xi = ±z f3 /2 (3.16) 

i=l 

where k is an arbitrary positive integer. We note that terms determined from the lower 
algebraic signs represent antipodal solutions. We also observe that although the optimal 
velocity increments must satisfy the form (I3.13p . and the sums on the left-hand sides must 
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match the specified end conditions on the right-hand sides, these velocity increments are 
otherwise arbitrary. In these singular cases there are non-unique optimal trajectories that 
satisfy the end conditions and minimize the related cost function (12.481) . We note that the 
left-hand side of ( 13 .161) is the minimum cost (12. 48ft . We note that for the singular solutions 
all velocity increments are stationary velocity impulses. In terms of the original variables 
r and 0, the singular solutions describe an outward (or inward) multi-impulse spiral from 
(0„, r„) to (9 f ,r f ). 

3.3 Three-Impulse Solutions 

For the remainder of this first paper we set k = 3 in (I2.58p - fl2.63p . 

According to (I2.6ip . for three- impulse solutions it is necessary that optimal velocity 
increments be applied where |p(0i)| — 1, i — 1,2,3. For primer vector equations (13. 61) and 
( 13. 7p this becomes 

A 2 [sin 2 0; + 4(cos0 4 - nf] = 1 (3.17) 

for i — 1, 2, 3. 

3.3.1 Stationary Solutions 

We first consider three-impulse solutions in which all the velocity increments are stationary 
velocity increments. These occur at values 9i where ( I2.62P and ( 13. 1 Tl) are satisfied. Using 
(13. 6 p and (13. 7p with ( I2.62p we find that if 0j is a stationary value then it is necessary that 

sin 9i cos 9i — 4 sin 0j(cos 0« — n) = 0. (3.18) 

This implies that either 1) sin0; = or 2) k = 3/4 cos 0j. 

We show now that the second alternative is not feasible because it violates (I2.63p . Using 
(13. 6p and (13. 7p and forming the function 

m = Pl (9f + p 2 (9) 2 - 1 (3.19) 

we find that if f'{9i) = where 0j satisfies 2) then /"(0j) > if cos 2 (0j) ^ 1. It so happens 
that cos 2 (0j) = 1 implies that the first alternative 1) is satisfied. Our argument shows that 
|p($)| > 1 in a neighborhood of 0j and equality holds in that neighborhood only at 0j if the 
second alternative 2) is satisfied. 

We therefore conclude that the first alternative 1) must be satisfied and 

9i = 7rz, i = 0,1,2,... . (3.20) 

Using two successive values of (I3.20p in (I3.17P we obtain 

2A(l-/t) = l (3.21) 
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2A(1 + k) = 1. (3.22) 
Clearly A = 1/2 and k = so the primer vector equations ( 13. 6p and ( 13.71) respectively become 

Pl (0) = ~sin0 (3.23) 

j9 2 (#) = cos^. (3.24) 

Forming the function (I3.19P we see that |p(0)| < 1 except at 9{ satisfying (I3.20p where 
equality holds, as seen in Figure 3. 

We shall pursue this problem further in our later examination of multi-impulse and de- 
generate solutions. It suffices to state at this point that if we substitute (I3.23P and (13.241) into 
(I2.60p . there are not unique solutions for a±, for k > 2. The three-impulse stationary 

solutions are therefore degenerate, and can be replaced by two-impulse solutions. 

In seeking non-degenerate three-impulse solutions, we must consider solutions in which 
all of the velocity increments are not applied at stationary values. 



3.3.2 Non-Degenerate Three-Impulse Solutions 

Throughout this paper we adopt the convention that the motion of increasing 9 is counter 
clockwise and that 9q < If ~~ 71 ^ < tt there are two cases, — 7r < 9 < and < 9 < tt. 
It is well known that there are at most four intersections of an ellipse such as (13. 5p with the 
unit circle. For the first case the intersections can be represented using at most four values 
9i, i = 1,2,3,4 where — n < 9i < n, i = 1,2,3,4. For the second case the four values 
representing intersections are considered on the interval < 9i < 2tt, % = 1, 2, 3, 4. 

We present the following lemma establishing necessary conditions for optimal non-degenerate 
three-impulse solutions. The primer vector loci for these solutions are depicted in Figure 1 
or Figure 2. The primer vector is defined by (I3.6p - (I3.8I) and the function / is defined by 

Lemma: It is necessary that an optimal non-degenerate three-impulse solution satisfy 
the following: 

9 1 = 9 , 9 3 = 9 f (3.25) 

82 = \(9 + 9 f ) (3.26) 

< 9 f - 9 < 2tt (3.27) 

f'(9 ) = -f(9 f ) < (3.28) 

f(9) < 0, 9 < 9 < e f . (3.29) 
There is an integer i such that either 2in < 9 < (2i + 1)tt and 

9 2 = (2z + 1)tt, and <p = 9 2 -(2i + 1)tt (3.30) 
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3 -3 
k=--(1-cos0 o ), — A < k < (3.31) 
8 4 

4 1 

A = - < A < 2 (3.32) 

5 + 3cos# 2 v ; 

or (2i — 1)tt <9q< 2iir and 

9 2 = 2in, and <fi = 9 2 — 2m (3.33) 

3 3 
K = -(l + cos0 o ), 0<k<- (3.34) 
8 4 

A = — ^, - < A < 2. (3.35) 

5 — 3 cos 9 2 

Proof: If a solution is non degenerate there must be a value 9 m on 9q < 9 < Of where 
6 m 7^ iii for any integer i and 9 m satisfies (I3.17p . otherwise all velocity increments are at 
stationary values and the solution is degenerate. For any integer j the value 9 n = 2nj — 9 m 
is distinct from 9 m and satisfies (13 . 1 7[) also. We let 9\ be the smaller of the numbers 9 m and 
9 n and 9 3 be the larger, and ^(^i + 9 3 ) = irj. 

There must also be a third velocity increment on a value 9 2 satisfying (13. 1 7)) where 
9q < 9 2 < 9f. For any integer I the number 2nl — 9 2 also satisfies (13.17|) . Since there can be 
only three impulses, this number is distinct from 9\ and ^3, and we must have 2nl — 9 2 = 9 2 , 
consequently 9 2 = irl for some integer I. 

We consider first the case where I is odd. We therefore write I = 2i + 1 for some integer 
i. Substituting 6\ and 9 2 into (13.171) and equating the left-hand sides, we obtain 

sin 2 9 1 + 4(cos 9 l - k) 2 = 4(1 + k) 2 . 
Expanding and simplifying we find that 

3(1 - cos 2 #1) + 8k(1 + cos 9i) = 0. 
This says that either cos^i = — 1 or 

3 

K = (1 — COS #1). 

8 

Since 9\ 7^ 7ri for any integer i, we have the latter. Differentiating f(9), evaluating the 
derivative at 9±, and substituting the above formula for k, we obtain 

= p(0i) ■ p'(0i) = -^A 2 sin^(l + cos^i). 

Next we observe that since 9\ + #3 = 2tij then 9% = 2nj — 9\ + <f> and 

/ / (03) = ^A 2 sin^ 1 (l + cos^ 1 ). 
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If sin^i < then f'(9x) > and f'(6s) < 0. This says that the primer vector exits the 
unit disk at 9\ and enters at 9 3 thus 9\, 9 2 and 9 3 cannot support a three- impulse solution 
since 61 < 9 3 and f !2.63|) must be satisfied. 

As a consequence of this argument, it is necessary that sin#i > hence f'(9i) < and 
f'(9 3 ) > 0. In order for a three-impulse solution to be supported on 9\, 62 and 9 3 and satisfy 
(I2.63P it is necessary that 9\ = 9 and 9 3 = 9f. It follows that sin# > resulting in the 
restriction 2m < 6 < (2i + l)n < 9f < (2i + 2)tt. (Recall that I — 2i + 1). The only integers 
that satisfy the restriction are I = j = 1, establishing (I3.25p - p.3ip for the first case recalling 
the definition of 0$. Substituting (I3.3ip into ( 13. 17ft and solving for A produces ( 13. 32 p . 

We consider next the case where / is even. We therefore write / = 2% for some integer i. 
Substituting 9\ and #2 into (I3.17P and equating the left-hand sides again, 



sin 2 9 X + 4(cos 9 1 - k) 2 = 4(1 - 



Expanding and simplifying as before we find 

3 

k = -(1 + COS^i). 
o 

Substituting this value of k into we get 

f(9 1 ) = p(9 1 ) ■ p'(0i) = ^A 2 sin^(l - cos^x) 

and similarly 

m)) = -^A 2 sin^(l-cos^). 

If sin^ > then > 0, f'(9 3 ) < and we again find that a three-impulse solution 

cannot be supported on 9i,9 2 and (9 3 and satisfy ( 12. 63ft because 9\ < 9 3 . 

It is therefore necessary that sin^i < 0, hence f'{9\) < and f'(9 3 ) > 0. In order 
that a three- impulse solution be supported on 61,62 and 6 3 and satisfy (12.63p . it is again 
necessary that 6\ = 6 and ^3 = 8f. This requires that sin^ < resulting in the restriction 
(2i — 1)tt < 9q < 2m < 6f < (2i + 1)tt. (Recall that I = 2i). The only integers satisfying 
this restriction are I — j — establishing ( I3.25p -( l3.29p . recalling again the definition of 6f, 
(I3.33P and (13.341) follow as well. Substituting (I3.34p into (I3.17P and solving for A reveals the 
formula (I3.35p . ■ 

Some of the statements in this theorem can be expressed differently by adding and sub- 
tracting 62 from 60 in ( 13.311) . ( 13.32p . (I3.34[) and ( 13.351) . We restate the lemma as the following: 

Corollary: It is necessary that a non-degenerate three-impulse solution satisfy (13.251) - 
(13T291) . and 

A = 4 (3.36) 

5-3cos(^) 
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Either 



3 

K = — — 



1 + cos( 



Of — , 



and the primer vector satisfies (13. 9 j) and (13. lOj) or 

3 



1 + cos( 



Of — O , 



(3.37) 



(3.38) 



and the primer satisfies (13.61) and (13. 7(1 . 

An immediate consequence of the lemma and the necessary and sufficient conditions 
(I2.58p -( l2.63p is the following fundamental theorem of non-degenerate three-impulse solutions. 

Theorem: The optimal impulsive rendezvous problem defined by ( 12.40(1 and (12.471) on the 
interval o < < Of has a non-degenerate three- impulse solution {Oi, 2 , 3 , AV 1; AV 2 , AV 3 } 
if and only if f l3T25l) -f l3T30l and ( 13T36H are satisfied, and either (ESD, (13TTU11 and <KWf\\ are 
satisfied or (13. 6ft . ( 13 . 7(1 (I3.33() and (I3.38() are satisfied, and 



AV ; 



-aip(0i), 



1,2,3 



(3.39) 



1=1 



cti > 0, 



O; 



1,2,3, 



(3.40) 
(3.41) 



where -R(#) is given from (12.521) and Zf from ( 12. 57(1 . 

From (13.311) . (13.32(1 . (13.34(1 and (I3.35P ) we see the bounds on A and k for which three- 
impulse solutions exist. These boundary primer vector loci are illustrated in Figure 3. 

We remark that an optimal solution in which (13.30l) -( )3.32jl are satisfied has an associated 
antipodal solution satisfying (13.33(1 - ( T3T35T) instead. Similarly if a generalized boundary point 
Zf has k satisfying (13.37(1 then its antipodal solution for the point— Zf has n satisfying (13.38(1 . 



3.3.3 Finding Non-Degenerate Three-Impulse Solutions 

We now apply the preceding theorem to find non-degenerate three impulse solutions. 

We consider the case where <p = 2 — 2iir and k and A are given respectively by ( 13.34(1 
and (13.351) . 0q and Of are known and 2 is known from (13.261) . Setting = — 2 the primer 
vector from (13.61) and (13.71) becomes 

p 1 (^) = Asin^ (3.42) 

p 2 (0) = 2\(cos0 - k) (3.43) 
where k and A respectively follow from (13 .34(1 and (13.351) with ^0 replacing 6?o- 
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Since (I2.49P is a fundamental matrix solution associated with A defined by (I2.47P so is 
also the matrix 



Analogous to (12.521) we define 



/ cos 6? sin (9 1 
— sin 9 cos 9 
1 



R{9) = ^-\9)B 



\ 



(3.44) 



and substitute into (I3.40p of the preceding theorem. Setting s 
0i = Xai, % = 1, 2, 3 we obtain 

±(5 -3s 2 -3c) |(5 -3c) ±(5 -3s 2 -3c) 



h(c-l) 



|(3 -5c) |(3c-5) 



ffl(l-c) 
5(3 -5c) 



01 \ 




( ~ z h 




H 




03 J 







(3.45) 

sin^ , c = cos^o an d 



(3.46) 



The region Zf is defined by the conical set generated by the column vectors of the 
coefficient matrix. To find the region in terms of inequalities, we first solve (I3.46p . use the 
identity c 2 = 1 — s 2 , and factor the denominators to obtain 



/ 




-1 



(5-3c)(l-c) 3(-s)(l-c) (5-3c)(l-c) 
-2(5c-3) g -2(3c 2 -3c+2) 



(5-3c) 2 (l-c) 



(5-3c) 2 (l-c) 




(3.47) 



-1 -1 -1 

\ (5-3c)(l-c) 3(-«)(l-c) (5-3c)(l-c) / 



It follows from (I3.27P and the definition of #0 that all of the denominators in this matrix 
are positive. Multiplying by the positive common denominators, employing ( 13.411) and the 
definitions of 0i, i — 1, 2, 3, we find 



3 (- s )^/i - ( 5 - 3c)^ /2 + 3(-s)zf 3 > 

(5c - 3)z fl + (3c 2 - 3c + 2)z h > 
3(-s)*/i + (5 - 3c)^ /2 - 3{s)z h > 



(3.48) 

(3.49) 
(3.50) 



From these inequalities it is seen that the subregion of Zf that admits three-impulse 
solutions described by fl3.33l) - fl3.35l) satisfies 



z h + z h > 



5 -3c 



\ z h 



(3.51) 
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(5c - 3)z fl + (3c 2 - 3c + 2)z fa > 0. (3.52) 

The antipodal region that admits three-impulse solutions by f l3.30p - fl3.32l) is found from 
the fact that <fi is decreased by ir from the preceding solutions, resulting in primer vectors 
having opposite signs. For this reason the left hand side of (I3.40p changes sign. The conse- 
quence is that all inequalities reverse in fl3.4Sp - fl3.52p . The two separate Zf regions admitting 
three-impulse solutions therefore have symmetry with respect to the origin. 

In the calculation of an optimal trajectory, the actual velocity increments Avj are cal- 
culated from the transformed velocity increments AVj. Solving for these from (12.371) and 
(E3HD we get 

Av 6i = ^AV 2i (3.53) 

i 2 

Av n = hiAVu + ^AV 2i (3.54) 

for i = 1, 2, 3. 

Example: Select 0„ = it/ 2 and 6 f = 3tt/2. From fl3T25|) and fl3T26|) we see that 6 X = vr/2, 
9 2 = vr and 6 3 = 3tt/2. We shall use fl3T33|) - (l3T35|) . If we use fl3T30|) - fl3T32|) instead we get 
a solution antipodal to this one. Using = tt — 2iix we obtain k = 3/8, A = 4/5 and 
§ Q = -tt/2. Eq. (JS3ZD becomes 

111 

01 = ~z h - -z h + -Zfy > (3.55) 

02 = ^z fl + ^z h > (3.56) 

3 = -z h + ~z h + -z h > (3.57) 

where (I3.5ip and f)3.52p become 

5 

z h + z h > ^\ z hl ~ 3z h + 2z h > °- ( 3 - 58 ) 

For the antipodal subregion the inequalities in (I3.58P are reversed. 
We observe that 

pi(0) = --sin0, p 2 (9) = --(cos9 + ~). (3.59) 
5 5 8 



Substituting into (13.391) we find 

AV, = n,A 

3/5 



AVx = ai ( A Jl ) (3.60) 



AV 2 = 




(3.61) 
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AV 3 = a 3 




(3.62) 



The antipodal solution reverses the directions of these velocity increments. The actual 
velocity increments follow from f !3.53|) and f)3.54p . The vector zj is calculated from f)2.50p 



These calculations were performed for the specific boundary conditions r = 8000km, 
v ro = —0.831929km/ sec, vq = 7 . 487362/cm / sec, andry = 6545.455/cm, v r , = —0.679267km/ sec, 
vg, = 7. 471940/cm/ sec. The result is 



which clearly satisfies (I3.58P and shows that the optimal rendezvous for these boundary 
conditions requires three impulses. The three velocity increments are found to be: Avg 1 = 
-0.23538km/ sec, and Av ri = -0.287688km/ sec at 9 1 = tt/2, Av 02 = -1.313191km/ sec 
and Av r2 = —7.691605 x 10~ 2 km/ sec at 9 2 = 7r where r(9 2 ) is approximately 7384/cm, and 
Ave :i = —0.193261km/ sec and Av r3 = 0.310108/cm/sec at 63 = 3tx/2. The fictitious nominal 
circular orbit has a radius of approximately 7272km. Simulation of the optimal three- 
impulse rendezvous trajectory connecting the initial and terminal conditions and showing 
the application of the three velocity impulses is depicted in Fig. 9. The maximum deviation 
of the rendezvous trajectory from the nominal circular orbit is approximately 10% of the 
nominal radius. 

4 Two-Impulse Solutions 

A two-impulse solution is a solution for which k = 2 in (I2.58p -f l2.63p . Recalling that / is 
defined by (I3.19P it follows that a non-degenerate two-impulse solution has exactly two zeros 
of / on the interval 9$ < 9 < 9f. 

The two-impulse solutions can be classified as either non-stationary solutions or station- 
ary solutions. 

The non-stationary solutions are non-degenerate and fall naturally into several categories. 
We recall that 9 = 9 + <fi where is an arbitrary constant. The following definitions apply in 
general to optimal k-impulse solutions but are especially relevant for k = 2. These definitions 
are illuminated by Figures 10 — 17. 

An optimal solution is called a non-stationary two-intersection solution if 9\ < 9 2 , 
f(9i) = f(9 2 ) = 0, f\9 1 ) ^ 0, f'(9 2 ) ^0,0<9 2 -9i< 2n, there are no other zeros of / on 
the interval 9q < 9 < 9f, and there is an integer i such that 9\ + 9 2 = 2iri. If < 9q < n 
then i — 1, if —n < 9 < then i = 0. Primer vector loci for two-intersection solutions are 
depicted in Figures 10 and 11. 



and fT2~57D . 



z f = (-2.777777 x 10~ 5 5.555555 x 10~ 5 ) T km 



-1 
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An optimal solution is called a non- stationary four-intersection solution if 6\ < 62, 

f(0i) — f{®2) = 0, < 9 2 — 9 1 < n, and there is an integer i such that m < 9% < 9 2 < + 
If — 71 < 6q < 7r then i = or i = — 1. Primer vector loci for four- intersection solutions can 
be seen in Figures 12 and 13. 

An optimal solution is called a non-stationary three-intersection solution if Q\ < 62, 
f(0i) = f(0 2 ) = 0, and either f'(0i) = or f'{6 2 ) = but not both. Examples of primer 
vector loci for three intersection solutions are presented in Figures 14 and 15. 

An optimal solution is called a stationary two-intersection solution if 9\ < 62, f{9\) = 
f(@2), f'(0i) = f'(@2) = 0, O2 — 0\ = vr and there are no other zeros of / on the interval 
0o < < Of- Primer vector loci of stationary two-intersection solutions are presented 
in Figure 16. If 82 — 6\ > it it is called a stationary multi-impulse solution and is 
degenerate. Primer vector loci of stationary multi-impulse solutions are depicted in Figure 
17. An optimal multi- impulse solution may have arbitrarily many impulses if Of — 6q is 
sufficiently large but always has an equivalent one-impulse or two-impulse solution having 
the same cost as will be shown in the subsequent investigation of degenerate solutions. 

The following theorem shows that there is a partition of the various categories of two- 
impulse solutions. 

Theorem: An optimal non-degenerate two-impulse solution is one and only one of the 
following: 

1. a stationary two- intersection solution 

2. a non- stationary two-intersection solution 

3. a non-stationary three-intersection solution 

4. a non-stationary four-intersection solution. 

Proof: Consider an optimal non-degenerate two-impulse solution having the primer 
vector given by (13. lip and (I3.12p . Forming the function ( 13.1 9p we obtain 

f{6) = A 2 [sin 2 9 + 4(cos 6 - k) 2 - 1. 

Differentiating this function, 

f(9) = -2A 2 sin£(3cos£-4fi;). 

If the velocity impulses are at 9\ and 62, then f{9\) = f{9%) = where 9$ < 9\ < 82 < Of. 
Since all solutions considered are two- impulse solutions then k = 2 in (I2.58p -f l2.63j) and there 
are no other zeroes of / on the interval 9 < 9 < 9f. 

1) If 9\ and 9 2 are integer multiples of n then f'(0i) = f'{9%) = and 9 2 — 0\ is an integer 
multiple of n. Since there are no other zeros on the interval 9q < 9 < 9f, then 92 — 0\ = n. 
This solution satisfies the definition of a stationary two-intersection solution. 
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2) If neither 9\ nor 82 are integer multiples of tt and there is an integer i such that 
0\ + 62 — 2/Ti then < 62 — 9\ < 2tt. If not, then 62 — 9\ > 2tt and there is an integer j 
such that either ff = 2ttj - 9 2 or & = 2ttj -9\. In either case 9 X < 9' < 9 2 and f{9') = 0, 
establishing more than two zeros of / on the interval 9q < 9 < 9f. 

We show now that f'(9i) 7^ and f'(9 2 ) 7^ 0. We argue by contradiction. Suppose 
f'{9i) = 0. Since 9 X is not an integer multiple of tt it is necessary that k = |cos0i. 
Differentiating f'{9) we find 

f"{9) = -2A 2 (6cos 2 #-4kcos#-3). 

Evaluating at 9\ and substituting for k into this expression we obtain 

fi(9 x ) = 6A 2 sin 2 ^. 

This is clearly positive since 9 X is not an integer multiple of tt. Since f{9i) = 0, then f(9) > 
on an interval where 9 > 9\. This violates (I2.63P establishing the contradiction. Similarly if 
we suppose f'{92) = we find that f(9)>0 on an interval where 9 < 9 2 , establishing the 
contradiction. As a consequence f'{9i) 7^ 0, /'(^2) 7^ 0, and the definition of a non-stationary 
two-intersection solution is satisfied. 

3) If 9i is not an integer multiple of 7r but (9 2 is an integer multiple of tt, clearly /'(^2) = 
and f'(9i) 7^ 0. If f'(9x) = the argument above establishes a contradiction. It follows that 
the definition of a non- stationary three-intersection solution is satisfied. 

If 9\ is an integer multiple of n but 9 2 is not an integer multiple of tt, clearly f'{9i) = 
and the above argument shows that f'(9 2 ) 7^ 0, so that in this case also, the definition of a 
non-stationary three-intersection solution is satisfied. 

4) Finally, if neither 9\ nor #2 are integer multiples of tt and if 9\ + #2 7^ 27ri for any 
integer i then there is an integer i such that iir < 9\ < ^2 < (i + 1)tt. The reason is as 
follows. If this statement is not true and if i is the largest integer such that in < 9\ then 
&2 > {i + l)vr. It follows then that there is an integer j such that if 9' = 2ixj — 9\ 7^ 9 2 or 
6' = 2nj — 9 2 7^ ^2 then 9% < 9' < ^ 2 and f(9') = 0. This establishes at least three zeros on 
9q < 9 < 9f and a contradiction. 

Having established that there is an integer i such that iir < 9\ < (9 2 < {i + l)ir we 
subtract 9\ in this inequality and obtain < (9 2 — 9\ < tt. The definition of a non-stationary 
four-intersection solution is therefore satisfied. ■ 

4.1 Non-Stationary Two-Intersection Solutions 

First we consider the category of solutions in which the primer vector loci demonstrate the 
geometry represented in Figures 10 and 11, the non-tangential intersection with the unit 
circle. In the following we specify that 9 < 9f. We define 9 r = 9 r — for any subscript r. 
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4.1.1 Fundamental Theorem of Non-Stationary Two-Intersection Solutions 

Theorem: The optimal impulsive rendezvous denned by (12.401) and (12.471) on the interval 
0q < & < Of satisfying the convention — ir < 0q < 7r has a non-stationary two-intersection 
solution {0i, 2 , AVi, AV 2 } if and only if 

e 1 = e , 6 2 = Of, (4.i) 

< f - O < 2vr, (4.2) 
/(0)<O, <0<0f, (4.3) 
f'(0o) = -f'(0f) < 0, (4.4) 



and either 



or 



and, in addition, 



4- 

< O < vr, 6 f = 2n- 0„, and = - J - vr, (4.5) 
- 7T < O < 0, f = -0 O , and = g ° + g/ , (4.6) 



AVx = -p(0i)ai, AV 2 = -p(0 2 )a 2 , (4.7) 
R{0 1 )p(0 1 )a 1 + i?(0 2 )p(0 2 )a 2 = -z f , (4.8) 
at > 0, a 2 > 0. (4.9) 
Proof: First we show that (14.ip - (l4.9p are necessary. 

If the optimization problem has a non-stationary two-intersection solution, then, < 
02 — 0\ < 27r where f{0i) = f{02) = and there is an integer i such that 0i + 2 = 2iri and 

< 0i < 02 < f . Since /'(0 a ) ^ and f(0 2 ) ^ 0, it follows from ([232]) that X = O 
and 2 = 0f and neither of these is an integer multiple of tt. From (I2.63P we have f(0) < 0, 
0o < < Of- Since there are not more than two zeros of / on this interval we must have 
f(0)<O,0 o <0<0f. We have established (144]) and (l4~3j) : (|4T4|) follows from the continuity 
of /, the fact that there are only two zeros of / on the interval 0o < < Of, f'(@o) 7^ 0, 
f'(0f) 7^ 0, and 0i and 2 = 2ni — 0\ can be substituted into the expression for f'(9). 

Since O is non- stationary we have two cases, either O<0 o <7ror— 7r<0 o <O. 

In the former case we observe that m = 2tt — o is a zero of / satisfying < O m — O < 2n. 
This requires that 0/ < m and (14 .2p is satisfied. Solving Of = 2n — 0q for cf) we find that 
(j) = ^ - establishing (I43|) . 

In the latter case we note that n = — 0o is a zero of / satisfying < 9 n — 9q < 2ir. This 
requires that 0/ < 0„ and (14. 2p is also satisfied in this case. Solving Of = — O for <f> we find 
that (f) = °^ 1 establishing ( 14. 6 p in this case. 

The expressions (14. 71) - (14. 9 p follow from the fact that an optimal solution satisfies (I2.58P - 
(I2.60P then ( 14 .3p shows the solution is non-degenerate and the inequalities in (14. 9 p are strict. 
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To show that (I4.ip - (l4.9p are also sufficient we note that they imply (58)-(63) of Reference 1 
with k = 2 showing optimality, where (14. 3 p shows that / has exactly two roots on 9q < 9 < 9f. 
The expressions (Hljl . flO]) and g3J show that + and f{9 2 ) + 0; gUJ and g2J or 

(14. 6 p show that 9\ + 9 2 = 2i\i where % = 1 or z = 0. 

These show that the optimal rendezvous is a non-stationary two-intersection solution 
completing the proof. ■ 

Remark: The convention — tt < 9 Q < tt is convenient as an assumption in this theorem 
but unnecessary. See the previous theorem for three-impulse solutions. Without this assump- 
tion (14. 5 p and (14. 6 p would assert the existence of an integer i such that 2iiT < 9q < (2i + 1)tt 
and (2i — 1)tt < 9 < 2m respectively. 



4.1.2 Finding Non-Stationary Two-Intersection Solutions 

Since A > we set f3\ = Acti, (3 2 = Aa 2 and (I4.9P becomes 

/?! > 0, (3 2 > 0. 



(4.10) 



We observe that is a fundamental matrix solution associated with A so R{9\) and 
R(9 2 ) m &y be inserted in (14. 8p . We select (14. 6 p so that 9f = —9 for the following develop- 
ment. There is no need to repeat this development for (I4.5P because it leads to the antipodal 
solution. 

Setting so = sin(9o) and Co = cos(9 ) (14. 8 p becomes 



61 

02 




( sl + 4c (c - k) Sq + 4c (c - k) \ 

-s c + 4s (c - k) s c - 4s (c - k) 

\ -4c (c - k) -4c (c - k) J 

For brevity we let 

g = 4(c - k) 

and obtain the following three equations in the unknowns 0i, (3 2 and go 



- z h 



So(c - go)(/3l - Aj) = Zfa 

go(/3i + /3 2 ) = z h . 

Solving these equations we find that 



1 



z h 



z h 



Sq + c g s {c - g . 



(4.11) 



(4.12) 

(4.13) 

(4.14) 
(4.15) 

(4.16) 
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h=-\ 



z h , z h 
_s 2 + c q s (c -q )_ 



(4.17) 



q = fe— . (4.18) 

Z /a + C Z /a 

These equations are subject to (I4.10p . The expression for k is found from (I4.12p and (I4.18p . 
Setting f(6 ) = we find 

\ = 2{Asl + ql) l l 2 . (4.19) 
Substituting fflOB"]) into ffl~T6l) and fT4TT7j) we obtain 



a = - ^L? * ( i + 1 ;° :\ ) ■ (4.2i) 



Z/ a + CqZ/3 / S Z/ 2 

2s o V coZfx + z h . 

In view of the inequality (" 14 .101) we find it necessary and sufficient for boundary conditions 
to satisfy 

z h + c z /3 < 0, |s Z/ 2 | < \c z fl + z h \ (4.22) 

to have non-stationary two- intersection solutions satisfying ( 14.61) . For antipodal solutions 
satisfying H4.5I) the necessary and sufficient conditions become 

z h +c z h > 0, \s z h \ < \c z h + z h \. (4.23) 



Example: Select 9 = — ir/2, Of = ir/2, c = and s = —1. From (14.181) we find 
that q = Zf z jzf x \ it follows from (I4.12p that k = — -J 3 -. We also find from (I4.19P that 

z 2 

A = 2(4 H — ir-) 1//2 although this formula will not be needed. 

For this example the region ( I4.22p of admissible solutions becomes 

z h < \ z h I < M- (4-24) 
The expressions (14.201) and (14.211) respectively become 

A = -f (l + |) (4.25) 

A = f 1 - (4-26) 



2 V 

The optimal velocity increments therefore follow from (14. 7p 

AVi = ( 1 ) A - Av * = ( _L ) * • < 4 - 27 » 
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The actual velocity increments are obtained from (13.531) and (I3.54p . For antipodal solutions 
the first inequality of (I4.24p and the directions of (I4.27P are reversed. 

This example was simulated for r = 6545km, v ro = —0.184km/ sec, 0q = 1.3xl0~ 3 rad/ sec, 
Tf = 8000A;m, v Tf = —0.710km/ sec, 0/ = 0.798 x 10~ 3 rad/ sec, boundary values which sat- 
isfy (23). The resulting optimal rendezvous orbit is depicted in Figure 18. The arrows in 
the figure indicate the application of the optimal velocity increments AVi = (AV x ,AV y ) = 
(0.4795, -0.8282) T km /sec, AV 2 = (AV X , AV y ) = (-0.0142, 0.2986) T km /sec. 

4.2 Non-Stationary Four-Intersection Solutions 

Next we consider the category of solutions in which the primer vector loci demonstrate the 
geometry represented in Figures 12 and 13, the four intersections. Again we specify that 
O < Of. 

4.2.1 Fundamental Theorem of Non-Stationary Four-Intersection Solutions 

The following theorem is fundamental for this type of two-impulse solutions. 

Theorem: The optimal impulsive rendezvous problem defined by (I2.40p and (I2.47P on 
the interval 0$ < < Of satisfying the convention —tt<0q<tt has a non-stationary 
four-intersection solution {0±, 2 , AVi, AV 2 } if and only if 

0i = 0o, 02 = 0/, (4.28) 

either < O < 0/ < tt or - tt < < f < 0, (4.29) 

and 

/(0)<O, O <0<0/, (4.30) 

f(0 o )<O, f(0/)>O, (4.31) 

3 3 3 

k = -(cos0 o + cos0/), --<«<-, (4.32) 
8 4 4 

A = [sin 2 O + 4(co8£o - ^f]- l/ \ \ < 2 (ll\ K \) ~ X ~ ( 3^ )V2 < 2 (433) 

and (H77)l - (0)l are satisfied. 

Proof: First we show that fl4T2l^ - fl4T33]) and ( H77|l -( IO]l are necessary. 

We suppose that the optimization problem has a non-stationary four-intersection solution. 
It follows that f(6 1 ) = f(6 2 ) = Q where / is defined by f l3~T9|) . O < 0i < 2 < f , there is 
an integer i such that in < 1 < 2 < (i + l)vr, < 2 - 1 < n, and k = 2 in (12. 58ft - ( 12331) 
from the definition of non-stationary four-intersection solution. By (12.621) there are no other 
zeros of / on O < < Of, hence 0i = 0o and 2 = 0f. From the convention — 7r < O < tt it 
follows that either i = and O<0o<0/<7rori = — 1 and — n < 0q < 0/ < tt. Since 0q 
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and 8f are the only zeros of / on 9q < 9 < 6f, (12.631) implies that /(#) < on the interval 
9 < 9 < Of. We have established fl2T28|) -( l2T30|) . 
Using fl3~Tll .f l3TT2l and ( 13191) we have 



f{9) = A 2 [sin 2 fl + 4(cos# 



K 



and 



f'{6) = -2A 2 sin^ (3 cos 9 -An). 



Setting f(9 ) = f(9 f ) we get 



(cos9f — cos6'o)[3(cos6 l o + cos9f) — 8k] = 0. 



Since < 9 f - 9 < vr then cos9 f ^ cos# establishing fT02|) . If f'(9 ) = then 
k = 3/4cos^o- Since 

f"(0) = -2A 2 (6cos 2 9 - 4k cos 9 - 3) 

then f"{9 ) = 6A 2 sin 2 9 > showing that f{9) > on an interval where 9 > 9q contradicting 
(I2.63P and establishing that f'{9 ) ^ 0. A similar argument shows that f'{9f) ^ 0. Either 
f{9 ) > or f'{9 f ) < implies that (14301) is false establishing ( 12~3T|) . 

The above arguments have shown (I4.28I) - (I4.32I) . Solving the equation f(9 ) = for A 
and minimizing and maximizing this expression with respect to 9q establishes (14.331) . Finally 
(HZD-dlSD follow from (12351) and (l2T5gj) - (l2^0j) by setting = 2. 

We now show that (I4.28l) - (l4.33p and (14. 71) - (14. 9 p are sufficient. These conditions show 
that (I2.58p - (l2.63p are satisfied for k = 2 resulting in an optimal two- impulse solution of the 
problem. This solution is non-degenerate because (14.311) implies that 9q and 9f are non- 
stationary. It is a four-intersection solution because (I4.28P and (I4.29P reveal that i = — 1 or 
i = and < 9 2 - 9\ < n follows from (14.291) after setting 9 X = 9 and 9 2 = 9f. ■ 



4.2.2 Finding Non-Stationary Four-Intersection Solutions 



Since A > we again set f3\ = Xaci, (3 2 = Xa 2 and (14.91) must be satisfied. For brevity we 
set Co = cos^O) sq = sin#o, Cf — cos#/ and Sf = sin9f. Since $(#) is a fundamental matrix 
solution associated with A we again use R(9i) and R(0 2 ) in (14.81) . Setting 9\ = 9$, 9 2 = 9f , 
substituting p(6 , ) and p(9f) into (14.81) and utilizing (I4.32p . after rearranging we obtain: 



( 1 + |c (c - Cf) 1 - |c/(c - Cf) \ 

\s Q {c Q -Cf) -\sf{c Q -Cf) 
\ -c - |(c - Cf) -Cf + Kco-Cf) J 




(4.34) 
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Solving these equations, we find that 



_ Sf 3(cq- C/ ) C f 
Pi - ~~£) h J) h 



S 3(00-0/) 

D Z * + D 



+ c 



if 



[9(g c + s f Cf) - 15(c s f + cjSq)} [9c 2 - 30c c f + 9c} + 16] 



6L> 



6L> 



where 



D = S + S f + -(C - C f ){c S f - CfS ). 



The velocity increments from ( 14. 7p become 

-s 



AVi 



-|co + \c t 



Pi, AV 2 



/3 2 



(4.35) 
(4.36) 

(4.37) 
(4.38) 

(4.39) 



In the case where < 6* < n, then D > and a non-degenerate four- intersection solution is 
found only in a Zf region where 



s f z h + 



3(c - c f ) 



c f 



z h < 0, s z fl - 



_3(c - c f ) 



+ c 



z h < 0, 



(4.40) 



and Zf 3 satisfies (I4.37p . Geometrically this says that there are solutions restricted to a sector 
described by (I4.40p of the plane described by (I4.37p . 

In the case where — n < 6q < then D < and the inequalities in (I4.40p are reversed but 
( I4.37P remains valid. Geometrically this reverses the sector in the plane with every vector 
Zf in the sector being replaced by — z f . This is the antipodal case. 

Of some interest is the special case where n = so that Co = — C/, so = Sf resulting 
in a primer locus having symmetry about a vertical line through the origin. The equations 
fl435|) - fl437j) simplify: 

z h z h 



01 

02 



Z h 



2(l + 3cg) 

z h 
2(1 + 3cl) 



4 

3s 



+ 



6s c 

z h 
6s co 



z h- 



(4.41) 
(4.42) 

(4.43) 



The velocity increments (I4.39P also simplify somewhat in this case. This type of solution 
exists only for the z/ region where 



3SqCq 3SqCq 

z fo < — 1 i o^2 z /i' z h > i , o^2 z h 



l + 3c^ 



l + 3c^ 



(4.44) 
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and Zf 3 satisfies (14. 43 p . consequently 



+ 

and 



^ n 3s ° c ° ^ ^ 3s ° c ° (A AK\ 



4s c 4s c , , 

2 /i < S Zf 3 < ~ - ; - j Z/i- (4.46) 



l + 3c 2^ — * - l + 3cg /1- 

Example: Select O = tt/6 and 0j = 57r/6, then cq = —Cf = a/3/2, so = sj = 1/2, 
A = -2/13z fl - 2y^/9z h , fa = -2/132/, + 2V3/9z h and 



i \ / i 



AVi = ( -V3 J ft - AV * ( A ) ft - (447) 

We select values of y an d Vf so that < and zj 2 = Zf 3 = 0. These conditions are satisfied 
by the boundary values r = 8000km, v ro = 0km/ sec, 0q = 8.37 x 10~ 4 rad/sec, T/ = 6545km, 
v Tj = —0.866km/ sec, Of = 1.251 x 10~ 3 rad/sec. Simulation of this example with these 
boundary values was performed, and the optimal rendezvous orbit is presented in Figure 19. 
The imposition of the velocity increments AVi = (AV x ,AV y ) = (0.0382, — 0.1216) T km /sec, 
AV 2 = (AV x ,AV y ) = (-0.0362, -0.1007) T km /sec is indicated by the arrows. 



4.3 Two-Impulse Non-Stationary Three-intersection Solutions 

Three-impulse non-stationary three-intersection solutions are discussed in Sec. 3.3.2. If Of 
is somewhat decreased or 0o increased from that of a non-degenerate three-impulse solution 
as shown in Figures 14 and 15 the result is a two-impulse non-stationary three-intersection 
solution. 



4.3.1 Fundamental Theorem of Two-impulse Non-Stationary Three-intersection 
Solutions 

Theorem: The optimal impulsive rendezvous problem defined by (I2.40p and (I2.47P on the 
interval 0q < < Of satisfying the convention — tt < 0o < tt has a two- impulse non- stationary 
three-intersection solution {0\, 2 , AVi, AV 2 } if and only if < 0$ < tt and 

9i = 0o, 02 = tt < f < 2tt - O , f(0o) <0,k = ^(cos0 o - 1), A = 4 (4.48) 

8 5 + 3 cos 0o 



or 



2tt - Of < O < X = tt, 2 = Of, f{0 f ) >0,k = %os0 f - 1), A = 4 (4.49) 

8 5 + 3cos#/ 
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or else — n < 9 < and 

9i = e J 2 = 0<9 f < -9 , f(6 ) < 0, k = ^(cos^o + 1), A = 4 , (4.50) 

8 5 — 3 cos O 

or 

_ 7T < < O < 0! = 0, 9 2 = 6 f , f(6 f ) > 0, k = \(caa6 f + 1), A = 4 , (4.51) 

O O — O COS Of 

and in either case 

f(e)<o,e <e<e f , (4.52) 

and (H77|) - fO|) are satisfied. 

Proof: First we show that (I4.48p - fl4.52p and (H~7]) - (H1?]) are necessary. 

If the optimization problem has a two-impulse non-stationary three-intersection solution 
then f{6 x ) = f(9 2 ) = where 9 < 6 l < 9 2 < 9 f and either f'(9 1 ) = or f(9 2 ) = 0a but not 
both, and — tt < 9q < it. We observe that 

f{9) = A 2 [sin 2 9 + 4(cos 9 - k) 2 ] - 1 

f{9) = -2A 2 sin 9(3 cos 9 - 4k) 
f"(9) = -2A 2 (6cos 2 fl -4k cos ^ -3) 
Equating f(9i) and 7(02) and performing some manipulations we obtain 

(cos 9 2 — cos 6*i) [3(cos 9\ + cos 9 2 ) — 8k] = 0. 

We shall show that cos#i ^ cos# 2 . 

Suppose that f'(9i) ^ and f'(9 2 ) = 0. If these are reversed the following argument is 
also valid with cyclic interchange of 9\ and 9 2 . 

Setting f'{9 2 ) = implies either sin^ 2 = or k = |cos# 2 but the latter implies f"{9 2 ) = 
6A 2 sin 2 9 2 > 0, which implies that f{9) > on an interval where 9 <9 2 because f(9 2 ) = 0, 
contradicting (I2.63p . It follows that sin# 2 = 0. If cos#i = cos# 2 then sin#i = also setting 
f{9\) — contrary to the original supposition consequently n = |(cos#i + cos# 2 ). 

Continuing f'{9 2 ) = implies sin# 2 — which implies that 9 2 = or 9 2 = it. 

If 9 2 = tt then < 9 X < 9 2 since f'{9 1 ) ^ 0. Setting k = 2 in fl232l) . we see that 
9i = 9 , so < 9 < n, 9 2 < 9 f , and 9 f < 2n - 9 because f(2n - 9 ) = 0. Clearly 
f'(9 Q ) < because f'(9 ) > with f(9o) = establishes an interval where 9 > 9 and (I2.63P 
is violated. Substituting cos^ 2 = —1 in the above expression for k reveals k = |(cos# — !)■ 
The expression for A comes from setting f(9 2 ) = substituting for k and solving for A. We 
have established f )4.48p . 

If 9 2 — the proof is analogous; — 7r < 9i < 9 2 since f{9\) ^ 0. Again setting k = 2 
in fr2T62|) we get 9 X = 9 so -tt < 9 < 0, 9 2 < 9 f and 9 f < -9 . Again f l2T63|) implies 
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that f'iPo) < 0. Substituting cos 6 2 = 1 into the above expression for k reveals that k = 
|(cos#o + 1)- Again setting f(9 2 ) = substituting for n and evaluating A completes (I4.50p . 

The expressions (I4.49P and (I4.5ip follow from f'(9i) = and f{9 2 ) ^ by repeating the 
preceding argument with 8i and 9 2 interchanged. These arguments lead to an interchange 
of 6q and Of and a reversal of all inequalities. 

The expression f l4T52|) is a restatement of fl2^3|) and g^-gl]) follow from f|238|) - (F2T60]) 
by setting k = 2. The inequalities in (14. 9p are strict because the solution is a non-degenerate 
two-impulse solution. 

To show that (Og|) -( 02|) and f^77jl - (jO]l are sufficient we note that they imply f l2"35j) - 
( I2.63P with k = 2 resulting in an optimal two- impulse solution. In each of fl4.48p - fl4.5ip 
h<h<0 2 < Of, and f(9x) = f(0 2 ) = and either f{0{) ^ or f'(0 2 ) ^ but not both, 
hence the solution is a non-stationary three-intersection solution and not degenerate. ■ 



4.3.2 Finding Two-Impulse Non-Stationary Three-Intersection Solutions 



First we consider the case where < Qq < tt, Q = Ox < Q 2 = tt < Of < 2tt — Q . Again we set 
Co = cos^o? s o — sin^o and observe that is a fundamental matrix solution associated 
with A. Substituting k = |(c — 1) into p(6 ) ) and p(^) the expression (14.81) becomes 



/ |(3 C 2 + 3c + 2) i(3c + 5) \ 



s (c + 1) 



\ -|(5c + 3) 
A solution exists if 



fa 



|(3c + 5) J 
2z f , 




(4.53) 



and 



3s (c + 1)' 
6g (c + l)z h - 2(3cg + 3c + 2)z h 
3s (c + l)(3c + 5) 

(3c + 5)z h 



3s 



(4.54) 
(4.55) 

(4.56) 



It is necessary from (14 .9p that (3i > and j3 2 > 0, and since the denominators in (I4.54p 
and (I4.55P must be positive, it is necessary that the Zf region satisfy (14 .561) and 



M+3co + 2)^ Zf >0 
h 3, (c + 1) ' Zh > U ' 



(4.57) 
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Geometrically, this region is a sector of a plane through the origin (56). If Zf is contained 
in this sector then (14 .7p determines optimal velocity increments 



AVi 



so 

5c +3 
4 



0i, AV 2 = - 





3cp+5 
4 



02 



(4.58) 



We note also that if — n < 9 < 0, 9 Q = 9i < 9 2 = < 9 f < —9o then we have an antipodal 
solution where zj is replaced by — Zf in fl4.53p - fl4.57p . 

Next we consider the case where 2tt — Of < < 0\ = it, 9 2 = Of- We set Cf = cos9f, 
Sf = sin9f and substitute k — |(c/ — 1) in p(#i) and p(#/). The expression (14. 8 j) becomes 



j ±(3c, + 5) \{3c) +3 C/ + 2) ^ 
§a/(c/ + l) 



/3i 
ft 



V|(3c / + 5) -|(5 C/ + 3) ; 
A solution exists if 




(4.59) 



01 



-6s f {c f + 1)^ + 2(34 + 3c f + 2)z h 







3s f (c f + l)(3c f + 5) 



3 5/ (c/ + l) 



and 



3c/ + 5 



(4.60) 
(4.61) 

(4.62) 



Since 0\ > 0, 2 > 0, and the denominators in f)4.60p and f)4.6ip are negative, The Zf region 
must satisfy (12.611) and 



z h < 



3c) + 3c f + 2 



z hi z h > 0- 



(4.63) 



3s/( C/ + l) 

Geometrically this is a sector of the plane fl4.62p . If Zf is contained in this sector, then 
optimal velocity increments are determined by (14. 7p 



AVi 





3c/+5 



U AV 2 = - 



5c/+3 



02 



(4.64) 



We note that if — it < 9 < 0, —Of < 9 < 9\ = 0, 9 2 = Of then we have an antipodal 
solution with Zf replaced by — zj in ( I4.59p - fl4.63p . 

Example: Select 9 = ir/2, 9 2 = it and Of = It follows that c = 0, s = 1, Cf = — | 
Sf = -^, (14.571) becomes 



z h > z h > 



(4.65) 
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and (I4.56P becomes 

5 

z h = z h ~ z z h- ( 4 - 66 ) 
Substituting into (57) of Reference 1 we obtain 

z h = -\vh + ~YVh + \vh + yo2, ( 4 - 67 ) 

v 7 ^ i Vs , Aco . 
z h = — 2~Vh ~ + ~Y y h - Vox + 2/o 3 , (4-68) 

z h = Vh-yvz- ( 4 - 69 ) 

Combining (ETM|) -( lZl)g|) . we find 

(2V3 - 6)y fl + (3^ + 2)y h + (-2y/3 + 3)y h + Ay 0l + 6y 02 - Ay 0s > 0, (4.70) 

- V3y fl - y h + V3y h - 2y 0l + 2y 0s > 0, (4.71) 
(3 - 5^3)^ - (3^3 + % /2 + (5V3 + 5)j/ /a - 10^ - 6y 02 + %, 3 = 0. (4.72) 

These inequalities are satisfied by the boundary values r = 7200km, v ro = 0km/ sec, 9 = 
1.033 x 10- 3 rad/sec, r f = 10375km, v rf = -3.hA7km/ sec, 6 f = 4.746 x 10- 4 rad/sec. This 
example was simulated for these boundary values, and the resulting optimal rendezvous 
orbit is presented in Figure 20. The optimal velocity increments, represented by the arrows 
on the figure are AV 1 = (AV X , AV y ) = (-0.3759, 0.5269) T km /sec, AV 2 = (AV x ,AV y ) = 
(-1.691, 0.9413) T km /sec. 



4.4 Non-Degenerate Stationary Two-Impulse Solutions 

A k-impulse solution is called stationary if f{9{) = 0, i = 1, . . . , k for distinct $i on the 
interval 8q < 9 < 9f implies f'{9i) = 0, i = 1, . . . , k. It is called non-degenerate if there is 
no /-impulse solution for I < k. It can be shown from the periodicity of (2.60) that stationary 
solutions for k > 2 are degenerate. 

Primer vector loci of stationary solutions are presented in Figures 16 and 17. Figure 16 
depicts a primer vector locus for a stationary non-degenerate two-impulse solution, whereas 
Figure 17 represents a primer vector locus of a stationary degenerate multi-impulse solution. 

Degenerate solutions are presented in a paper that follows. 



4.4.1 Fundamental Theorem of Non-Degenerate Stationary Two-Impulse So- 
lutions 

We shall accept the convention — 7r < 9q < tt. If #o is outside this interval, a simple change 
of variable can be employed. 
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Theorem: The optimal impulsive rendezvous defined by ( I2.40p and (I2.47P on the interval 
$o < 9 < Qf satisfying the convention — 7r < 9 < tt has a non-degenerate stationary two- 
impulse solution {6*1, # 2 , AVi, AV 2 } if and only if 

k = 0,\ = ± (4.73) 
n<9 f -9 < 3tt, (4.74) 

and if — n < 9q < then 

0! = 0, 9 2 = tv, <p = 9 1 = 9 2 -n, tt < Of < 2vr, (4.75) 
but if < 9 < 7r then 

0i = tt, 9 2 = 2n, <p = 9 X - tt = 2 - 2tt, 2tt < /; (4.76) 

and in either case 

f(Q) = -- sin 2 9 < 0, 9 <9<9 f (4.77) 

and (j£7j) - (0|) are satisfied. 

Proof: First we show that (I4.73jl - p.77p and f l4.7p -( )4T9l) are necessary. It is shown in 
(13 . 191) that a stationary solution must be supported on integer multiples of 7r; non-degeneracy 
requires this support at the two smallest multiples of tt greater than or equal 9 , consequently 
9f must be placed after the second smallest multiple of it and previous to the third. The 
results dHZU)- < K7Q\f follow and fT4\73|) is found from fl3T20|) and < K21\f . Substituting fl473|) 
into the expression for the primer vector, we have 

f{9) = ^ sin 2 9 + cos 2 9 - 1 

which results in fETTTD . We obtain fl4"?7] -( l4T9"]l by setting k = 2 in (1238D - (123(1 : since the 
solution is non-degenerate the inequalities in (14. 9p must be strict. 

It is seen that (I4.73p - (l4.77p and (I4.7p - (l4.9p are sufficient because they show that this 
solution is stationary and f !2.58[) - (I2.63j) are satisfied for k = 2. This solution is non-degenerate 
because the inequalities in (14. 9p are strict. ■ 

We now apply this theorem in order to find non-degenerate stationary solutions. We 
investigate two cases, — n < 9 < and < 9 < n . 
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4.4.2 Finding Non-Degenerate Stationary Solutions 



Case 1: Suppose — 7r < 6q < then 9\ = 0, 62 = vr and n < 9f < 2tt. Observing that 
p(^ 1 ) T = (0, 1) and p{6 2 ) T = (0, -1), using R(9) from f[232l) . we see that gSD becomes 



/ 2 2\ 




a: 2 




(4.78) 



V-2 2 J 

The solution of this system of equations is 



«i = —^ z h ~ z h 



«2 



where 



From (14. 9p . it is necessary that 



z h = 0. 



(4.79) 
(4.80) 

(4.81) 
(4.82) 



z h < ~\ z h\- 

Optimal two-impulse solutions for this case are found for boundary conditions satisfying 
(14.81 j) and (I4.82p . Geometrically the z/ region is restricted to a sector (I4.82p of the plane 
(I4.8ip . If these boundary conditions are satisfied, the optimal velocity increments are applied 
at 9 1 = 0, 82 = vr and are calculated from (14. 7p : 



AVi 







AV, 





a 2 



(4.83) 



Case 2: Suppose < 6 Q < ir then B\ = 71,62 = 27r and 2tt < 9f < 3tt. We note that 
p(^) r = (0, -1) and p{6 2 ) T = (0, 1). In this case dHHJ) becomes 



/2 2 \ 




en 1 

«2 




(4.84) 



V 2 -2/ 
The solution of this system of equations is 



"1 = ~^i z h+ z h) 



(4.85) 
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a 2 = -~(z fl -z f3 ) (4.86) 

where again (14.81 j) is satisfied. From (I4.9p . it is necessary also for this case that (I4.82p is 
satisfied. 

Optimal two-impulse solutions for this case also follow from the boundary conditions 
satisfying (I4.8ip and (I4.82p . For these boundary conditions the optimal velocity increments 
are applied at Q\ — ir, and 2 = 2ir and are calculated from (14. 7p : 

A Vl = ( ° ) , AV 2 = ( ) . (4.87) 

Example: Select 9o = 0, 9f = n so that Q\ = /9 , 9 2 = 9f and = 9 = 9f — ir. Using 
( 12301) and fl236l) with 9 in ( 14301) we have 

Z h = ~Vfi + Vh ~ VOi + VOa 

z h = -Vh - vo 2 = o 

The condition ( I4.8ip becomes 

Vo 2 + Vh = 

and (14.821) becomes 

-Vh + Vh - yoi + yo 3 < -\vf 3 — 2/o 3 1 - 

If yf 3 > yo 3 this becomes 
but if yf 3 < yo 3 it becomes 

2yo 3 < yoi +y/i, 

and if yj 3 = yo 3 then 

z h = %o 3 - (yox +yh) < o. 
The boundary values r = 7200/cm, v ro = Okm/sec, 6 = 1.15 58 x 10~ 3 rad/sec, rf = 7200km, 
v Tf = Okm/sec, Of = 1.1554 x 10~ 3 rad/ sec. satisfy the first of the last three inequalities, and 
are used for simulation of this example. The resulting optimal rendezvous is shown in Figure 
21. The optimal velocity increments AVi = {AV X , AV y ) = (1.7756 x 1(T 3 , -0.8779) T 'km/ sec, 
AV 2 = (AV x ,AV y ) = (-4.3162 x 1(T 4 , -0.8038) T fcm/sec are indicated by arrows in the 
figure. 



5 One-Impulse Solutions 

For certain boundary conditions a minimizing solution may contain only one impulse. This 
situation is investigated here. 



41 



5.1 Fundamental Theorem of Non-Degenerate One-Impulse Solu- 
tions 

We apply (I2.58p -( )2.63p where k = 1. If a one-impulse minimizing solution {0i,AVi} is 
non-degenerate then AVi ^ 0. For this case applications are facilitated without use of the 
vector Zf. 

Theorem: The optimal impulsive rendezvous problem defined by (I2.40p and (I2.47P has 
a non-degenerate one-impulse solution {0i, AVi} if and only if 

0x = O or f'ipx) = or 9 X = 0/, (5.1) 

f(Q x ) = and f(9) < on O < 9 < 6 f , (5.2) 

ori = |AVi|, (5.3) 
AVi 

-pie,), (5.4) 



|AVi| 

fiAVx = $(0x - 9 f )y f - $(0i - 0o)yo- (5.5) 

Proof: Setting k = 1 and AVi ^ in (EMJ, flZHjgp and fl23ll -f l2T63l . these expressions 
become equivalent to fl5.1l) - fl5.4p . In f)2.60p we replace 0! by 0i obtaining _R(0 1 )AV 1 = zj. 
Multiplying on the left by $(0i) in view of (12.521) the expression (I5.5P emerges. ■ 

5.2 Finding Non-Degenerate One-Impulse Solutions 

The expression (15.11) is resolved into two situations, the non- stationary solutions where Q\ = 
0o or 0i = 0/, or the stationary solutions where f'{9\) = 0. 

5.2.1 Non-Stationary Solutions 

In the following we set 6 d = Of — O , c& = cos0d, and Sd = sm0 d . 
Case 1: Suppose 0! = O then (15.51) becomes 

£?AVi = ${-6 d )y f - y . (5.6) 
Utilizing (12.401) and (I2.4ip and solving (15.61) we obtain 

AV n = SdVh - °dVh + Vo 2 , (5-7) 

AV 12 = ~(y h -y 03 ), (5.8) 

where 

Vo 1 = c d y h - s d y h + y h . (5.9) 
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Expressions (I5.7p - (l5.8p are necessary and sufficient to have one-impulse solutions supported 
at 6*0. 

Case 2: Suppose 0\ = Of then (I5.5P becomes 

BAVi = y/ -$(^)y . (5.10) 

Solving, 

AVLi = -y /a - s d y 0l + c d y 02} (5.11) 

Ay 12 = -i(y /3 -i/03), (5.12) 

where 

Vh = CdVot + s d yo 2 + Vo 3 - ( 5 - 13 ) 

Expressions (I5.11I) - (I5.13I) are necessary and sufficient to have one-impulse solutions supported 
at f . 

5.2.2 Stationary Solutions 

In the following we set c = cos#o> s o — sin# , Cf = cosOf, and Sf = sin9f. 
Case 1: Suppose 0\ = and — n < 9 < < Of < it, then (15.51) becomes 

BAY, = *{-6f)y f - *{-e )y . (5.14) 

Solving, noting ( 15. 4ft and the fact that = if 0\ is a stationary point, we obtain 

AV n = 0, (5.15) 

&V 12 = -^(y h -y 03 ), (5.16) 

where 

Co2/0! - s y 02 + y 0s = c f y fl - s f y h + y h , (5.17) 

SoVoi + c y 02 = s f y fl + c f y h . (5.18) 

Expressions (15. 1 5 f) - ( T5TT8T) are necessary and sufficient to have stationary one- impulse solu- 
tions at 0\ = 0. 

Case 2: Suppose 0\ = tt and < < n < Of < 2n, then (I5.5P becomes 

BAV, = $(tt - f )y f - $(tt - o )y o . (5.19) 

Solving, as before we find 

AV U = 0, (5.20) 

AV 12 = ~(y h -y 03 ), (5.21) 
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where 

CoZ/Oi - So2/o 2 ~ Voa = c fVh - s fVf2 - Ufa ( 5 - 22 ) 
SoVo 1 - c y 02 = s f y h - c f y h . (5.23) 

Expressions fl5.20p - fl5.23p are necessary and sufficient to have stationary one- impulse solu- 
tions at §1 = Ti. 

Conclusions 

A planar impulsive rendezvous problem can be solved for initial and terminal positions and 
velocities near those associated with a nominal circular orbit. The rendezvous trajectories 
produced are exact for this restricted two-body problem and approach optimality as the 
boundary conditions approach boundary conditions of a nominal circular orbit. 

By replacing the cost function f)2.1ip by the related cost function f)2.48p . it was found 
that velocity impulses that minimize 02.48)) produce a close approximation of a minimum of 
(12. lip if the boundary conditions are near those of a nominal circular orbit. 

The boundary conditions are incorporated into a vector Zf defined by (12.571) which de- 
termines the number of velocity impulses required and from which these velocity increments 
can be calculated. It was found that a non-degenerate rendezvous trajectory requires three 
impulses if and only if the boundary conditions are such that z/ or —Zf satisfies f)3.5ip and 
f )3.52p . It was found that the proper placement of the velocity impulses is at the initial, termi- 
nal, and mid-point values of the true anomaly. The calculations are simple enough that the 
three velocity increments could be found from a hand calculator that contains trigonometric 
functions, if necessary. An example of a three-impulse rendezvous problem with boundary 
conditions that lead to quick solution was presented. Simulation of the rendezvous trajectory 
that resulted showed that throughout the flight the deviation of the radial distance remained 
within ten percent of the radius of a nominal circular orbit. 

It was found that non-degenerate two-impulse solutions of an optimal impulsive ren- 
dezvous near a circular orbit can be classified according to the ways the locus of a primer 
vector can intersect the unit circle. These optimal two-impulse solutions fall into four cate- 
gories. 

For each of these categories necessary and sufficient conditions for an optimal solution 
were presented. Distinct regions of the boundary conditions that admit optimal two-impulse 
solutions for each category were displayed. A closed-form solution of the optimal velocity 
increments for each category can be found. Some examples and simulations were presented. 

Necessary and sufficient conditions for optimal non-degenerate one- impulse solutions were 
also found. These optimal one-impulse solutions consist of four types. Distinct regions of 



44 



boundary conditions admitting the optimal one-impulse solutions were also displayed. For 
each type a closed-form optimal velocity increment can be found. 

Although this work emphasizes non-degenerate rendezvous in the vicinity of a circular 
orbit, it also presents a framework for studies of stationary, degenerate, singular rendezvous 
and optimal rendezvous beyond the vicinity of a circular orbit. Presentation of results of 
these studies is to follow. In these results we will show that degenerate and singular solutions 
are more than interesting curiosities. Not only do they provide important understanding to 
the subject of impulsive rendezvous, but they also comprise important solutions to some 
impulsive rendezvous problems. 
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List of Captions 



Fig. 1 - Primer locus and unit circle for three- impulse solution. 

Fig. 2 - Primer locus and unit circle for antipodal three-impulse solution. 

Fig. 3 - Primer locus and unit circle for a type of two-impulse solution 

Fig. 4 - Primer locus and unit circle for another type of two-impulse solution 

Fig. 5 - Primer locus and unit circle for antipodal two-impulse solution 

associated with Fig 4. 
Fig. 6 - Primer locus and unit circle for a type of one-impulse solution. 
Fig. 7 - Primer locus and unit circle for another type of one-impulse solution. 
Fig. 8 - Primer locus and unit circle for yet another type of one-impulse solution. 
Fig. 9 - A three-impulse rendezvous trajectory. 

Fig. 10 - Primer vector loci for two-intersection solutions where < O < n, + 0f = 2n. 
Fig. 11 - Primer vector loci for two-intersection solutions where — n < O < 0, 0o + 0/ = 0. 
Fig. 12 -Primer vector loci for four-intersection solutions where O<0 o <0/<tt- 
Fig. 13 -Primer vector loci for four-intersection solutions where — n < O < Of < 0. 
Fig. 14 - Primer vector locus of a two-impulse three-intersection solution where 
< O = 0! < 7T = 2 < B f < 2?r - o . 

Fig. 15 - Primer vector locus of a two-impulse three-intersection solution where 
< 2vr - 9 f < O < 7T = 0i < 2 = Of- 

Fig. 16 - Primer vector locus for a stationary non-degenerate two-impulse solution where 

-71 < O < 0, 7T < Of < 27T, 71 < f - O < 3tt. 

Fig. 17 - Primer vector locus for a stationary degenerate multi-impulse solution having 
an equivalent two- impulse representation where < 0o < 0/ or — 7r<0 o <0/, 
3tt < f - O . 

Fig. 18 - Rendezvous orbit (red) of a non- stationary two-intersection solution 
O = -7t/2, Of = tt/2, r = 6545A;m, v ro = -0.784km/ sec, O = 1.3 x I0~ 3 rad/sec 
T f = 8000km, v rf = -0.710km/sec, 0/ = 0.798 x 10~ 3 rad/sec 

Orbits satisfying initial conditions (green) and terminal conditions(blue) are displayed. 
Fig. 19 - Rendezvous orbit (red) of a four-intersection solution 

O = tt/6, Of = 5tt/6, r = 8000A;m, v ro = 0km/ sec, O = 8.37 x lO^rad/sec 
Tf = 6545/cm, v rf = -0.8QQkm/sec, f = 1.25 x lO^rad/ sec. 

Orbits satisfying initial conditions (green) and terminal conditions(blue) are displayed. 

Fig. 20 - Rendezvous orbit (red) of a two impulse nonstationary three-intersection 
solution O = tt/2, 2 = ix, Of = 47r/3, r = 7200km, v ro = 0km/ sec, O = 1.033x lO^rad/ sec 
Tf = 10375km, v rf = —3.547km/ sec, Of = 4.746 x 10~ 3 rad/ sec. 

Orbits satisfying initial conditions (green) and terminal conditions(blue) are displayed. 

Fig. 21 - Rendezvous orbit (red) of a nondegenerate stationary solution 
O = 0, Of = 7r, r = 7200km, v ro = 0km/ sec, O = 1.1558 x lO^rad/ sec 
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Tj = 7200km, v rf = Okm/ 'sec, Of = 1.1556 x 10 3 rad/sec. 

Orbits satisfying initial conditions (green) and terminal conditions(blue) are displayed. 
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Rendezvous orbit of a non-stationary two-intersection solution 
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Two Impulse Three Intersection Rendezvous 
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Figure 20: 
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